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Abstract 



We present full renormalisation group equations for the MSSM with R-parity viola- 
tion, including all soft supersymmetry breaking terms. The inclusion of dimensionless 
R-parity violating couplings can generate many possible low energy effects through 
the generation of off-diagonal soft masses violating lepton and quark flavour, and 
through the generation of lepton-Higgs mixing. We discuss the relation between the 
weak and unification scale R-parity violation and study the effects on neutrino mass 
generation and — > 67. 



1 Introduction 



Among the different proposals to extend the Standard Model (SM), supersymmetric 
models are probably the most attractive ones as their fundamental ingredient, su- 
persymmetry (SUSY) |jl|, provides a very elegant solution to the so-called hierarchy 
problem. However the most popular SUSY extension of the SM is not the most general 
one, as it is commonly assumed that a discrete symmetry, known as R-parity, forbids 
all the baryon {B) and lepton (L) number violating couplings, hence avoiding possible 
rapid proton decay. To each field is assigned an R-parity given by R = (_ i)3-B+l+25 
{S being the spin) so that all the SM particles have R = +1 while their SUSY part- 
ners have R = —1. Immediate consequences of imposing such a symmetry are that 
SUSY particles are produced in pairs and that the lightest SUSY particle (LSP) is 
stable, therefore becoming a dark matter candidate. 

However, there is no obvious reason why B or L violating interactions should not 
be allowed separately to break R-parity @, |^ , and the apparent constraints through 
the requirement that any new physics should not ruin baryogenesis in the early uni- 
verse have been shown to be very weak |Q. The corresponding terms in the superpo- 
tential are: 



where u, d, e stand for the u-type, d-type and lepton singlet superfields, and Q, L are 
the quark and lepton doublet superfields respectively, which generate dimensionless 
couplings in the Lagrangian, given by: 



= 2^ijkUidjdk 



(1.1) 



and 



Wji^ — -XijkLiLjCk + X[ji^LiQjdk , 



(1.2) 




(1.3) 
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and: 



+ Kjk {'^LidLjd.Rk - eRiULjd*Rk + J^LidLjdRk - eiiULjdRk 
+ VLidLjdRk - eLiULjdRk) + h.c. 



;i-4) 



Besides these, we also expect dimensionful terms analogous to the R-parity conserving 
soft breaking terms; that is, we get soft trilinear couplings and masses which violate 
B or L : 

soft = l^'ijk^idjdk + h-c. + J2 Kb'Pa'fb (1-5) 

a,b 

Vii^soft = ^CijkLiLjek + C[jj,LiQjdk + h.c. + ^ ml^(pa^b + DiLiH2 + h.c. (1.6) 

a,b 

Here the soft masses shown include such terms as ^^i.^i which violates L, and also 
terms violating lepton and quark flavour symmetries (where each generation is as- 
signed its own Li and Bi) whose effects will be in addition to the contribution from 
the CKM sector. For a discussion of our notation, see Appendix A. 

R-parity violation has been considered by many authors in the past and recently 
there has been an increasing interest in the renormalisation group properties of R- 
parity violating couplings 0, §|, 0; if- This is specially interesting when we try to 
embed the SM as the low energy limit of a more fundamental Grand Unified The- 
ory (GUT). In such a framework we have a prediction for the values of the different 
couplings at the unification scale, Mqut, and we can therefore use their renormal- 
isation group equations (RGEs) to obtain the corresponding values at low energies 
and compare them with the existing limits. This applies to both dimensionless and 
dimensionful couplings, for which the RGEs are given (the latter for the first time) 
in Appendix A. 
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Concerning the limits which can be imposed on these couphngs, there have been 
many attempts to constrain them through their impact on rare processes in the lab- 
oratory [y, |TU], |rT| and it turns out that flavour changing neutral current (FCNC) 



processes in particular provide us with a very powerful tool to test physics beyond 
the SM. The presence of extra particles, SUSY partners in this case, together with 
lepton and quark flavour violating couplings, including here off-diagonal soft masses, 
in our fundamental theory can induce transitions forbidden in the pure SM and, 
moreover, strongly contrained by experiment. Consequently an upper bound can be 
derived on different products of dimensionless couplings and even on some of the di- 
mensionful ones. Furthermore, some of these newly generated couplings can give rise 
to neutrino masses and mixings which are again experimentally constrained. 

In general we expect that once we allow the violation of lepton and quark flavour 
symmetries through any operator, including the R-parity violating dimensionless cou- 
plings shown in equations (|1.2|) and ( |1.3|) , other such operators will be generated 
through the RGEs. In particular, we here study for the first time how the R-parity 
violating couplings generate soft terms which can lead to large FCNC effects in the 
lepton sector process /z — >■ 67. 

Therefore, the aim of this paper is to present a detailed study of the limits which 
can be obtained on these R-parity violating couplings at the GUT scale by imposing 
FCNC and neutrino mass bounds at low energy. In particular, we have focussed on 
L violating interactions, but we plan to study B violation in a forthcoming paper. 

In section 2 we discuss the dimensionless R-parity violating couplings. We sketch 
analytical solutions for the case when the gauge contribution dominate in their RGEs 
and obtain triviality bounds for different values of the strong coupling 03. In sec- 
tion 3 we clarify which basis is the most appropriate to deal with lepton and quark 
flavour violation and also the lepton-Higgs (LH) mixing which arises from R-parity 



3 



violation. Section 4 is devoted to a study of sneutrino VEVs, their generation from 
the dimensionful RGEs, and the hmits which can be set on them by using the present 
constraints on neutrino masses. In section 5 we put hmits on various R-parity vi- 
olating couplings at Mqut by running them to the electroweak scale and requiring 
that the full SUSY contribution to the FCNC process /i — > 67 be less than its actual 
experimental bound. We perform a full numerical study and give some analytical 
approximations. Finally in section 6 we present our conclusions. 

2 Dimensionless Couplings 

Before turning to the dimensionful couplings, it will be helpful to discuss what conclu- 
sions we can draw from the RGEs in Appendix A with respect to the dimensionless 
ones. We start by remarking that there have recently been a number of analyses 
@) S % § studying the RGEs for some or all of the dimensionless couplings. Our 
RGEs are consistent with both references and 0, apart from a couple of ty- 
pographical errors in the gauge contributions to the RGEs in reference 0. More 
complete lists of RGEs are given in reference [0 , but here we disagree on a number of 
coefficients, even at one loop. Our equations give a different sign for the third term 
in equation (2.13) and the last term in (2.16) of that paper. Finally, we agree with 
the equations of reference p. 

While the equations are too complex to be usefully solved analytically in general, 
we might expect that often many (or all) of the R-parity violating Yukawa couplings 
will be tiny relative to the gauge couplings, and hence that only the effect of the latter 
will be relevant. In this limit, it is straightforward to solve the RGEs analytically ]12[ 
as follows. 

aa(to) 



aait) 



l-26,(t-to)«a(to) 



\{t) = (l-22(t-to)ai(to))'' (l-2(t-to)«2(to)) 
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x(l + 6(t-to)«3(to)) 
\'{t) = (l-22(t-to)«i(to))''^''(l-2(t-to)«2(to))''''' (2.1) 



x(l + 6(t-to)a3(to))" 



4/6 



A"(t) = (l-22(t-to)«i(to))'"^''(l-2(t-to)a2(to))''^' 



:(l + 6(t-to)a3(ito 



476 



Here we have used the convention that t = \nfi where fi is the MS renormahsation 
scale, while 

ha = (11, 1, -3) 

Ca = (3, 3, 0) (2.2) 
cl = (7/9, 3, 16/3) 
c'i = (4/3, 0, 8) 

and the label a runs over gauge groups 1,2,3. Although it is possible to derive an- 
alytical solutions even when the Yukawa couplings are large [|12|, in practice the 



expressions become so complicated that they are not particularly useful. 

If we solve these equations with t = Iii^Mgut/Mz), and using Mqut = 2 x 
lO^^GeV, as = 0.11,0.12,0.13, we find 

X{Mz) = 1.5 \{Mgut) 

X'{Mz) = 3.4, 3.6, 3.7 \'{Mgut) (2.3) 
X"{Mz) = 4.0, 4.3, 4.7 X"{Mgut) 

These results break down at the 10% level only when A(M^), A'(M^), \"{Mz) exceed 
0.27, 0.15, 0.15, which are very much larger than the values which we shall consider 
in our FCNC and neutrino calculations. For the R-parity violating couplings which 
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couple directly to ht (and, for sufficiently large tan/3, to hb) the ratio of the low energy 
to high energy value is reduced somewhat, typically by up to around 20%. 

Another interesting feature, and one which has been discussed in a number of 
recent papers f^, 0, §, is that it is possible to derive triviality bounds on the R- 
parity violating couplings in exactly the same way as for the top quark Yukawa. Such 
bounds depend quite strongly on the values of tan/3 and a^. We give bounds derived 
from our RGEs in Table 1 by simply setting the coupling to be very large at the GUT 
scale, either with small ht or ht also at triviality, and running down to a scale of order 
the top mass. These results agree quite well with those of references [^, || , and agree 
also with those given in reference ^ where applicable. 



Coupling 


Limit, as = 0.11 


Limit, as = 0.12 


Limit, 0^3 = 0.13 


A 


0.92 


0.92 


0.92 


A' 


1.08 


1.10 


1.12 


A" 


1.17 


1.20 


1.23 


A 


0.92 (1.08) 


0.92 (1.11) 


0.92 (1.13) 


y 


1.08 (1.08) 


1.10 (1.11) 


1.12 (1.13) 


^x3x 


1.00 (1.00) 


1.02 (1.03) 


1.04 (1.04) 


^xxx 


1.17 (1.08) 


1.20 (1.11) 


1.23 (1.13) 


X" 


1.02 (0.93) 


1.05 (0.95) 


1.08 (0.96) 



Table 1. Triviality limits on various couplings derived firstly by considering ht small, and secondly 
in the case where ht approaches its triviality limit. In the latter case we show the low energy value 
of ht in brackets. Here "x" can represent either 1 or 2 (or 3 if the index corresponds to an L or e 
and we are not in the very large tan p regime) , and the calculation has been performed for a 

variety of different as^Mz). 



3 Choice of Basis 



We shall be interested in deriving constraints on R-parity violation which arise 
through lepton flavour violation (LFV), in which processes are allowed to violate 
lepton number generation by generation, quark flavour violation (QFV), or the mix- 
ing between leptonic and higgs superflelds (LH mixing). Given the nature of these 
effects, it will be convenient here to discuss our choice of basis, and how this basis 
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affects our calculations. Indeed, it has recently been pointed out that bounds derived 
in a quark mass eigenstate basis can be rather different from those derived with the 
same couplings in a weak eigenstate basis |T^. While the physics is invariant, it is 
important that we choose a convenient (and consistent) basis to work in. 

It is well known that, by judicious rotations of the fermion fields in the Standard 
Model, it is possible to work in a basis where the mass eigenstates are as far as 
possible weak eigenstates, thus ehminating FCNC at tree level and ensuring lepton 
flavour conservation at all orders even through charged current processes. In terms of 
supersymmetric parameters, the three fermion mass matrices correspond to Yukawa 
matrices h'^j, h^j, hfj of which the first two can be simultaneously diagonalised. Given 
the usual assumptions of universality at the GUT scale, the soft mass matrices are 
diagonal while the trilinear terms are proportional to the corresponding Yukawas, 
and so there is no LFV, and the only QFV effects beyond the usual SM ones are 



those coming from the CKM matrix, which are relatively small [|14[ |T^. This picture 
can break down when non-universal soft masses are allowed at the GUT scale, for 
example because one is considering a realistic unified theory, but this is outside the 
scope of our analysis. 

It is clear from the RGBs presented in Appendix A that for some arbitrary choice 
of R-parity violating couplings there will be LFV, QFV, and LH mixing effects. The 
last case, for example, is driven through diagrams of the type shown in Figure 1. 
These generate non-diagonal anomalous dimension matrix elements of form 7^.^ which 
violate lepton number, and analagous diagrams exist which conserve L and B but 
generate LFV or QFV. With the insertion of appropriate spurion terms for the soft 
masses we will also find LFV and QFV effects in the soft sector. 

Thus, even working in a basis where the Yukawa matrices which generate the 
fermion masses are diagonal at the GUT scale, this will not necessarily be true at the 
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electroweak scale. It is simplest then to invoke the GIM mechanism and to perform 
a field rotation so as to eliminate the most awkward effects and simplify calculations, 
which we find to be the quark and lepton mass eigenstate basis. The question now 
is what effect this rotation will have on the LFV, QFV and LH mixing terms in the 
soft sector. 

Continuing with the case of LFV, in general this field rotation must take the form 
L — > [/L, e — > V*e, suppressing all indices for simplicity, so that 

Lh^'eHi L{U^h^V*)eHi (3.1) 

with U and V unitary matrices selected so that (U'^h^V*) is diagonal. We now 
consider the effect of this rotation on the slepton mass matrices. We find that 

ml ml := U^mlU (3.2) 

If we assume that 

m| = X 7 + 5m^ (3.3) 

where / is the identity matrix and 6m^ is a matrix with the only non-zero elements 
small and off-diagonal, while similarly U — l + 0{e) with e parametrising the rotation 
and assumed to be small, we can see that 

m| = mi + 0(e5m'^) (3.4) 

Clearly, if we have sufficiently small off-diagonal elements in the mass and Yukawa 
matrices, we are quite safe in taking the non-diagonal terms in the mass-matrix to be 
the same before and after we have rotated back into the appropriate mass eigenstate 
basis. It is straightforward to check that this argument is valid for all the LFV effects 
which we consider. 

The argument above assumes degenerate masses for the relevant generations and 
so will clearly be inaccurate in the case of any process involving the third generation, 
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such as 6 ^ 57, and for the case where we consider LH mixing. Here there is no 
alternative but to perform the rotation into the appropriate basis exphcitly. The case 
for LH mixing is somewhat different, since here the preferred basis for simphcity is not 
a mass eigenstate basis but one in which there are no couphngs in the superpotential 
of form fiiLiH2, but the same principle applies. As will be shown explicitly in the 
next section, even continuously rotating the fields so as to ensure that /Xj remains zero 
at all scales will not remove all the LH mixing effects in, for example, the running of 
the Yukawa couplings, the soft bilinear terms, and Tn\.^^. 

We conclude by noting that in the fermion mass eigenstate basis the non-diagonal 
soft masses will appear in the form of flavour violating mass insertions in our diagrams, 
rather than in the vertices themselves. 



4 Sneutrino VEVs 
4.1 General Discussion 

One signature of R-parity violation which has been much discussed is that of sneutrino 
VEVs. These typically arise because of the existence of /ij, Dj, and fn\.Hi terms which 
explicitly cause the effective potential to contain terms linear in the sneutrino field, 
either from explicit P, [l^ or spontaneous fl^ R-parity violation (the origin of R- 



parity violation is irrelevant if it appears in this sector, and recently the implications 
of such effects have been considered in the context of GUTs ||18[). They can also 
be caused by one loop effects involving dimensionless R-parity violating couplings 



m 20 



Once sneutrinos have acquired VEVs, neutrinos and neutralinos mixQ, so that we 
may derive bounds on R-parity violating terms by imposing experimental limits on 
neutrino masses. In fact, it has recently been shown that it is possible not merely 



^As in fact, charginos and Icptons, and Higgscs and sneutrinos 
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to bound R-parity violating terms, but to select such terms so as to explain the 
present rather complicated array of experimental measurements of neutrino masses 
and oscillations |^ . 

The potential which we must minimise takes the form 

+2{m]j^L, + fX4l^i)uik + 2Diku2 (4.1) 

Here li is the VEV of Lj, with z, j running from 1 to 3, we drop all terms of order 
and above (since we expect I <^ u). We define as usual z/j =< Hi >, with u"^ = 
uf + z/| = (174GeV)2 and tan/3 = z/2/z/i. 

In order to study these effects, it is most convenient to follow reference and 
perform a field rotation mixing the Li and Hi such that only one of the /ij, which 
we choose to be /i4, is non-zero. Doing so will not merely change the dimensionful 
couplings, but will also generate new contributions to the dimensionless couplings A 
and A' from hr and ht. In this basis, however, both the Higgs-sneutrino potential and 
the neutrino-neutralino mass matrix take on relatively simple forms. Up to order l/u, 
and assuming that the off-diagonal terms in the mass matrix Tn\.i^_ are small, we find 
that 

' mi^ + iM|cos2/3 
There are two dangerous approximations which we have made and which should be 

justified. The first is that we have only considered the Higgs potential at tree level. 

This is well known not to be reliable in general, but we would not expect that the 

large radiative corrections due to quark and squark loops should greatly affect the 

calculation of sneutrino VEVs. Hence we expect our analysis to be reliable so long 

as we derive the Higgs VEVs z/i, z/2 from the full one- loop effective potential before 

substituting them into equation ( ^4.2| ). 
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Our second approximation is that we have not considered the processes discussed 
in detail in reference in which one-loop diagrams generate sneutrino VEVs even 
with all dimensionful R-parity violating couplings zero. We shall ignore this point, 
and simply draw the reader's attention to the fact that the effects we study will in 
general be supplementary to these. 

The neutralino-neutrino mass matrix then takes the form, in the basis of 2-spinors 
given by [iB^, iWl U], 



I 



Ml 






M2 



TT^^i'^i TS^i'^a -775^1^* \ 



72^?2^^1 






.V2i 



72^2^^2 -^92 














(4.3) 



suggesting that the neutrino mass generated is of order {g\ + gl)lf/2M, where M is 
some typical neutralino mass. 

The simplest constraints on neutrino masses for the first generation are the di- 



rect bounds as quoted for example in reference ||22|. We could also use experimental 
bounds on neutrinoless double beta decay, on neutrino oscillation, and from cosmol- 
ogy. The oscillation effects are rather complicated and we shall not consider them 
in this analysis except to note that the latter would give constraints on two or more 
R-party violating couplings considered together, similarly to the case where the neu- 
trino masses and mixings are generated directly from the dimensionless couplings, and 
that in any case products of As are already very tightly constrained by these effects 
|19| , pO| , p3| . In light of the complexity of our results, we shall not derive bounds as 
such, but instead will discover how large the typical impact of the R-parity violating 
couplings can be. 
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4.2 Generation of Sneutrino VEVs from the RGEs 



Given the discussion of sneutrino VEVs above, we wish to investigate whether the 
RGEs can tell us anything about the likely size of such effects. The usual scenario has 
supersymmetry breaking in the hidden sector communicated to the visible sector only 
by gravity, so that the soft breaking terms take the form of trilinear and bilinear terms 
proportional to the corresponding Yukawa couplings and /i terms, with universal soft 
masses. If we make this assumption, then at the GUT scale (or perhaps more plausibly 
the Planck scale) we can simply rotate away the //j for i 7^ 4 as described above, and 
in doing so we will also set all = except for D4 without generating any non-zero 
rri^j^i^.. Hence we shall assume that at the GUT scale we have done this, and so our 
initial configuration will have R-parity violation only through dimensionless couplings 
and trilinear terms. 

If we then run our initial conditions down to the weak scale, then we will find that 
here sneutrino VEVs will only be generated through the following non-zero products 
of couplings. 

\izi,K or A-gg/ife generating /i^, A, ^H^Li 

Ci33hr or C-gg/ife generating Di (4.4) 

Ci33VT or Cl^^rjb generating 171%^^ 

These of course are exactly what we would expect, since in order to generate sneutrino 
VEVs we must have terms mixing Lj and Hi, as shown in Figure 1. 

Once we run down to the weak scale, however, we must perform another change 
of basis to that in which = 0, and this will introduce additional effects. The 
simplest way to describe these is by recasting the relevant RGEs in the basis where 
we continuously rotate as we run so as to ensure that 

In this basis we find that for i 7^ 4 and assuming that the diagonal terms in the soft 
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mass matrices are very much greater than the off-diagonal terms and //4 ^ //j 

167r2— ~ -Xi3shr{2ml + 2ml + 2ml^) (4.6) 

-2Cj33r?r - 6C-33?76 

Note that although it is possible to cancel off some of the L^-Hi mixing terms by 
a field redefinition, not all terms can be cancelled simultaneously. Here for example 

the rotation has removed the terms from the RGE for 7^1/^^,., but increased the 
coefficient of the mi,, terms. 

Wc sec that the effects are likely to be largest when tan/3 is large, giving large 
hb and /i^, but that the dependence is clearly rather complicated. More significantly, 
the sneutrino VEV will in general be proportional to the R-parity violating coupling, 
and hence the neutrino mass to the coupling squared. 

4.3 Numerical Results 

Many of the numerical methods and assumptions which we use here will be standard 
for calculations throughout this paper. We have as input parameters rrit, a^, mo, Aq, 
Ml/2, tan/3 = 1^2 /I'l and the sign of /i. We run the Yukawa and gauge couplings up to 
the GUT scale, where we impose universality of the soft masses, but not unification, 
and set any R-parity violating couplings which we wish to investigate. ttIq and M1/2 
are the values of all the diagonal elements of the scalar soft mass-squared matrices and 
the universal gaugino mass respectively. We set each trilinear coupling equal to the 
corresponding Yukawa coupling multiplied by Aq, so that for example C[jj^{MGUT) — 
X[jj^{MGUT)Ao- The masses and couplings are then run down to low energy to give 
output. and D4 are chosen so as to give the correct minimum of the one-loop 



13 



effective potential. We do not include any threshold corrections at either scale, using 
the non-SUSY RGBs below M1/2 and the SUSY RGBs above it. 

Typical results are shown in Figure 2a, where we set Xi33{Mgut) = 0.01 and Figure 
2b where instead we set XissiMcuT) = 0.001. Input parameters are rrit = 175GeV, 
asi^Mz) = 0.12, Aq = 0, M1/2 = SOOGeV, and we plot the absolute magnitude of the 
sneutrino VEV and the corresponding neutrino mass as a function of m^^, the soft 
mass of the left-handed slepton, for both signs of fi and for tan P =2,20. We note that 
in certain cases the sneutrino VEV changes sign at some value where the contributions 
for the non-zero Di and ^^iLi exactly cancel. It is clear from equation (O) that. 



for given Aq and for one or other sign of /i4, Di and Tn'\j^i^^ will have opposite sign, 
and here we see that complete cancellation can occur when > 0. Here larger tan /3 
gives a larger sneutrino VEV and hence larger neutrino mass. 



Although we might naively expect from the form of equation (^4.3| ) that increasing 
Ml/2 (for fixed w^lJ should cause the resulting neutrino mass to decrease, increasing 
Ml/2 also changes the value of ^4 and affects the running of all the dimensionful 
parameters and so the situation is rather complicated. In Figure 3 we plot the absolute 
magnitude of the sneutrino VEV and the resulting neutrino mass against M1/2 for 
the same parameters as Figure 2b, but with rriLi held fixed at SOOGeV. Again we see 
that the relative cancellation occurs to give zero neutrino mass for some set of input 
parameters, and is noticable that even with M1/2 as large as SOOGeV, the neutrino 
mass is still increasing with increasing M1/2 in at least some cases. 

The effect of varying Aq is also large, as shown in Figure 4. The main limit on 
the size of Aq is that its absolute value should be less than around 3 times a typical 



slepton mass to avoid charge or colour breaking minima p4|. In Figure 4 we vary 
Aq over the whole of this range for the same parameters as in Figure 2b, but with 
rriLi = Ml/2 = SOOGeV. From this figure it is clear that the cancellation of the effect 
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can occur for either sign of /i, depending on the signs and relative magnitudes of the 
low energy value of the trilinear couplings. 

We now turn to a discussion of how these calculations can be used to derive a 
limit on the couplings. It is clearly straightforward to derive bounds on Ajsa and X[^^ 
for a point on one or other of Figures 2a to 4 by calculating which value of A or A' 
gives a neutrino mass which is acceptable in the light of the experimental limits, given 
that the sneutrino VEV is directly proportional to A133 (or A'^33), and so the neutrino 
mass generated is proportional to the square of each of the coupling. What is perhaps 
most interesting here is that although the constraint tends to become less tight with 
increasing values of the soft masses, as is generally the case for the constraints on 
R-parity violating couplings, the decrease can be slow relative to the case for limits 
derived from one-loop diagrams with internal sfermion lines, as these are almost 
always power-law suppressed by soft masses. We typically find that even for soft 
masses of order SOOGeV the bounds are relatively tight, with XissiMcur) ~ 0.01 and 
X[^^{Mgut) ~ 0.001 both giving neutrino masses of order IkeV or more, suggesting 
a bound on these couplings which is at least an order of magnitude tighter for the 
first generation. These constraints are much tighter than those in the literature on 
individual couplings [|, |Tl|, since normally the tightest constraints are on products 



of couplings. It should however be noted that we are quoting limits on couplings at 
the GUT scale, which should be converted to values at the weak scale using equation 
( p.3| ) for comparison with those given by most other authors. 

Such limits are obviously greatly hampered by the dependence on the many input 
parameters, on the assumptions about the GUT scale structure, and in particular by 
the fact that there are strong cancellations between terms. We note that one possible 
cause of partial cancellation, namely having A^sa and A^gg of the same order and 
opposite sign, is unlikely in the context of certain GUTs, where we would typically 
expect Xiss^McuT) = Ksai^Gur) just as hb^McuT) = hriMcur)- Hence our results 
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should be regarded as indicative of the hkely magnitude of the effect rather than being 
able to strongly constrain the electroweak scale values of the couplings in isolation. 
Nevertheless, it is significant that the RGEs typically give such large effects, and 
in constructing any model which includes any of the couplings which can generate 
neutrino masses it is necessary either to set them extremely small or else to fine-tune 
away the unwanted effects. 



5 /i ^ 67 

5.1 Contributions 

The process /i — 67 is one of the most tightly constrained examples of FCNC Since 
this process cannot occur in the SM, its observation would be pressing evidence for 
new physics beyond the SM. Experiments have calculated an upper bound on 
the branching ratio (BR) of 

BR{fi 67) < 4.9 X 10"^^ at 90% CL (5.1) 

In SUSY models, a non-zero rate can be generated through non-diagonal slepton mass 



matrices [|I4|, |23] , and also through the direct effects of R-parity violating couplings 



0, |16|, [T^. However, as noted above, R-parity violation induces LFV through soft 
terms, and so in any realistic model where the latter effects occur, so will the former, 
and here we shall consider the two effects together. 

The total branching ratio for the process — >■ 67 can be written as: 

BRifi ^ e7) = ^ {\Alr\' + \AnL\') (5.2) 

where Gp is the Fermi constant, and Alr, Arl are the total amplitudes to the LR 
and RL transitions respectively. Here 

Ai = Af + if + if'" (5.3) 
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with i — LR, RL and the different terms are as follows. The various contributions to 
^Li?) ^LR arc shown in Figures 5 to 7 respectively. Diagrams contributing 

to are not shown, but are similar to the LR case for and the neutralino con- 
tributions to A'^i^. Both A^]^ and the chargino contributions to A^"^ are proportional 
to the electron mass and thus neglected. 

Results for the various amplitudes are : 
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sin^ df.^ cos^ 9e 
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(5.4) 



where e is the electromagnetic constant, is the mass of the corresponding sneu- 
trino, m-,(i),(2) denote the two mass eigenvalues of the mass matrix for the gen- 
eration of sleptons and sin^e^.; cos^e^ a-re elements of the corresponding orthogonal 
diagonalisation matrix. These are given by: 
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SO that state e*^^) is defined by 



(5.5) 



e*-^'' >= cos^ekz, > -|-sin6'e|eK > 



(5.6) 



The second amplitude is given by: 



A^' 

^LR 



+ ^^^(F.(.5?)+2F.(.<?))-^(2F,(4'') + f.(4'')) (5.7) 
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^ '(1).(2) 
5 -^ij 



J(l).(2) 



The various functions are mostly found in and |2^, and are reproduced in Ap- 
pendix B for convenience, while x^P'*'^'' 

Finally the third amplitude is: 
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where sw (cvk) is the sine (cosine) of the Weinberg angle, and g is the SU(2) gauge 
coupling constant. Also, 
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and the mixing matrices f/, l^, A^, A^' are defined as in ||29|. Finally is the 



trilinear coupling for the second generation of sleptons. We use Am? = Am? 



LfJ-L 



and Am? - instead of m? r and m^ ^ as in Appendix A to avoid confusion between 
generation and mass eigenstate labels. 

The RL amplitudes are given by: 
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with A^i^ = and: 
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In calculating these amplitudes, we have adopted the usual approximation that the 
electron mass is negligible relative to the muon mass, and we work to first order in 
both m^/Mw and m^{A^ + /itan/3)/(m|^ — m|^). 

Finally we note that we have neglected many diagrams which are suppressed by the 
mixing between neutralinos and neutrinos, or between charginos and charged leptons 
1^, IB]. Two typical examples of such diagrams contributing to A^ji are shown in 



Figure 8. Such terms will exist due to mixing of neutral fermions as shown by the 
mass matrix of equation (|4.3|) , and similarly from the mass matrix for charginos and 
sleptons which contributes to the lagrangian (in two-spinor notation) 

/ M2 g2iy2 



[-iW~,H^ ,eu 




(5.13) 



921^1 /W4 hik 
\ Qik hiVi 

Here we have presented only one generation for simplicity, and hi is its mass- generating 
Yukawa coupling. However, even for the third generation the sneutrino VEV U must 
be quite small to avoid an excessive neutrino mass, and hence the mixing is small. 
Relative to the direct proportional contributions shown in Figure 5 such diagrams 
with a predominantly lepton line are suppressed by glU/i^XM) (where a factor §2/^ is 
from the differing couplings, and g2li/M is from the mixing angle). Thus we expect 
this diagram to contribute less than the direct contribution except when A is so small 
that both are negligible, in which case we would anyway not expect the generation of 
large /j, since from our earlier discussion of sneutrino VEVs we always have li < XM 
by at least an order of magnitude. Similarly, if the internal fermion line is predomi- 
nantly neutralino or chargino, such contributions are similar to those shown in Figure 
7 but suppressed by a relative factor /jA/M, and so negligible. 
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Further diagrams can be drawn with mixing between sneutrinos and neutral Higgs 
states, and between charged sleptons and charged Higgs states, but by similar argu- 
ments we are justified in neglecting those too. We note however that although we 
expect such diagrams to give weaker bounds that for the case where corresponding 
diagrams without sneutrino VEV insertions exist, they can provide very weak con- 
straints on products of couplings which we cannot otherwise bound. For example 
Figure 8a will allow a bound on the product Ai33Ai2i. 

5.2 Analytical Discussion 

Given the complexity if the expressions shown above, we can try to see whether in 
certain limits we are able to derive analytical bounds on the R-parity violating cou- 
plings. In each case, we shall set only two R-parity violating dimensionless couplings 
non-zero at Mqut and see what effects they generate. We begin by discussing the 
relatively simple "direct" contributions to A'^ and A'^' before turning to the "indirect" 
contributions of A^*". 

Due to the antisymmetry of the Xijk couplings with respect to the first two indices, 
we will always have i = 3 in equation (|5.4|) and therefore m^. = m^-^. Also note 
that, for the first two generations, the weak interaction and the mass eigenstates 
are practically the same, so we can take then sin^ 6'ej. = 0, cos^ 9e^. = 1 and also 
"^^(1) (2) = fn'^-L R, with k = 1,2. This gives us a bound on products of As at the 
low-energy scale 

\X3ik{Mz)X32k{Mz)\ < 4.6 X 10-^ (yOOG^J k = l,2, (5.14) 

assuming that m^^ ~ mg^ ~ m. Bounds for the case k = 3 are pretty similar 
to these, though their dependence on the mixing in the stau sector makes them less 
straightforward. However in the case of maximal mixing (that is, sin^ 9^ = cos^ Or = \, 
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and with one stau very much hghter than the other) we can see that: 



|A3i3(Mz)A323(Mz)| <2.3 X 10- 



/ lOOGev y /lOOGevV 



(5.15) 



where rrif^ is the hghter stau mass eigenvalue. Therefore we see that there is a 
possibihty of canceUation inside the bracket, thus considerably reducing the branching 
ration generated by this process and hence relaxing the bound. A further problem 
with these analytical limits on slepton mediated diagrams is that we are assuming 
that the right and left handed sleptons have the same masses, and in practice this 
assumption can break down completely, as will become apparent when we come to 
consider our exact numerical results. 

We turn now to possible limits on these couplings coming from the RL contri- 
butions to /i — s> 67 given in eq. (^.4]). In this case we have one pair Ajji, Xij2, with 
i,j = 1, 2, and i j which satisfy the bound 

\k,i{Mz)X^j2{Mz)\ < 2.3 X 10-^ (yOOG^) ''^ ^ ^^'^^^ 



Note that there is a relative factor of two between equations (|5.14| ) and ( |5.16|) from 



the fact that we must set Xiji = —Xjn and so cannot set only one A 7^ at once. The 
remaining cases, i = 3, j = 1, 2, allow cancellations in a similar way to the case of 



equation (|5.15|). 



Concerning the A' couplings, we can turn back to eq. ( [5.7] ) and notice that for 
i,j = 1,2 again the weak interaction and mass eigenstates are the same, so that 



sin^6'„^ = sin^6'd^ = 0, cos^ 6'„^ = cos'^ 9d^ = 1 and "^|(i),(2) = ^Il,r, "^j(i),(2) 



^jL.fl. Furthermore, due to the small values of the quark masses for these first two 
generations, Xji and x[j are very small as well, and Fi{x), F2{x) can be replaced by 
1/6, 1/12 respectively. All this allows us to deduce that 

\X'u,{Mz)X'2.,{Mz)\ < 4.6 X 10"^ (yOOG^) = ^'^ ^^'^^^ 
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where here we assume that the squark masses are degenerate at m. However when 
either i, j or both are equal to 3, everything is going to depend on the mixing between 
stops and/or sbottoms; in that case we find a tighter bound (a factor of 2 smaller that 
for the other couplings) when that mixing is large by assuming that the contribution 
of the lightest mass eigenvalue (sbottom for i = 1,2 and j = 3, and stop for the rest) 
will be the dominant one in the corresponding amplitude. 

So far we have presented analytical bounds on products of couplings at the elec- 
troweak scale. However, note that these can immediately be translated into bounds 
on the same couplings at Mquti just by performing the corresponding running up 
from Mz- Moreover, the approximate solutions to the RGEs sketched in eq.( p.3|) will 
provide the reader with automatic bounds for most of the relevant cases. 

By looking at eqs. ( ^.8|) , ( ^.12|) it is clear that not much can be extracted from 



them analytically except in limits so simple as not to be useful. However, we can at 
least remark upon their sign and order of magnitude. If we consider the corresponding 
RGEs for LFV with diagonal mass terms dominant, 

mniC mnj^ (5.18) 



-\-2Ci.fYifiCjn 



mn 

+^CimnC^jmn) ) (5.19) 



we can solve these equations only in the limit where we only consider the contributions 
from the largest gauge coupling after imposing unification at Mqut- However, it is 
more useful simply to give the empirical results below : 

ml,^ {Mz) = - E >^mniXmnj{MGUT){l.7ml + CGM^^ - 0.3AoAfi/2 + 0.6v42) 

mn 
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mn 

+ E ALnA;>„„(MG^^)(15m^ + 30Mi% - 2OA0M1/2 + 5A^) (5.20) 

mn 

The accuracy of these formulae is roughly at the 10% level for the case where we deal 
with A only, but is much lower for the case with A' due to the strong dependence on 
the QCD coupling, whose value is uncertain and which typically appears in the RGEs 
with larger coefficients. In the latter case we find that the M^yg coefficient can vary 
by almost a factor of two, and so the results should be regarded as merely indicative 
of the likely scale of the terms generated. 



Hence we see that the terms Am? - / m? and Am? ,- / m? are of the same order 
as XiijX2ij + 'iX'njX'2ij and XijiXij2 respectively where the couplings are taken at Mz, 
and are rather larger when we have non-zero A's than As. Given the range of values 
of the various functions and the other couplings, it is then straightforward to check 
that we expect the neutralino and chargino contributions can be at least comparable 
in magnitude to the direct contributions, and often larger, particularly for the case 
of the A' where the internal particles in the loop for the direct case are heavier but 
the impact from Am^ is larger than for the A case. 

We also note that except H{x) > G{x), L{x) > F{x), and hence if all the mixing 
matrix elements are significant we expect the second term in each of the chargino 
and neutralino contributions to dominate. Thus we might expect that usually A^"^ 
will grow with increasing tan/5, and fall off only slowly with increasing chargino and 
neutralino masses, since the explicit mass factor will partially cancel the dependency 
in the function. 

Another important point which we should make here is that when we add up all 
the contributions to Alr from equations ( ^.4|) , (ISTTI), and (|5.8|) , and to Arl from 



and ( p.l2| ), to obtain the total, there are likely to be cancellations. For example, if we 
consider XiijXi2j > in the limit where there is very little gaugino-higgsino mixing, 
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then the contributions to and the first chargino term of are both positive, 
while the contribution to the first neutrahno contribution to A'^^ is negative, and 
other terms will vanish. 

5.3 Numerical Results 

In any case, we can derive through these equations new bounds on products of A 
and A'. As discussed above for the sneutrino VEV calculation, we simply impose 
unification at Mqut with some pair of R-parity couplings non-zero, evolve down to 
low energy and then investigate the amplitudes for yU 67. Given that here there 
will be sensitivity to virtually every parameter through the complicated dependence 
on the spectrum, we shall attempt to pick typical scenarios and examine whether 
reasonable conclusions can be drawn. 

Firstly, let us consider the impact of X'ui^Mgut) = ^^iA^gut) = 0.001. Here we 
note that we expect A^33 to be negligible because of its generically very large impact on 
sneutrino VEVs and hence neutrino masses as discussed above, while we shall avoid 
A'^j3 and A'^3j for simplicity. For mo = M1/2 =100GeV, we expect that the values 
chosen will give a contribution to the amplitude which is less than the experimental 
limit by about two orders of magnitude using Eq. ( ^.17|) . 

We show the resulting contributions to the amplitude as a function of M1/2 in 
Figure 9. Here we have set tan/? = 10, mg = lOOGeV, Aq = 0, /i4 > 0, and 
\'iii{Mgut) = '^'2ii{^gut) = 0.001. What is remarkable about this figure is that it 
is clear that in this case the direct contributions are completely negligible relative to 
those from the chargino and neutralino mediated diagrams. Similarly keeping M1/2 = 
lOOGeV, so that the chargino is always light, and varying mo, we find contributions of 
form of Figure 10. Here we find a curious cancellation between neutralino diagrams, 
so that the neutralino contribution at one point changes sign relative to the chargino 
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contribution, but is in any case always small. Finally we show the dependence on 
in Figure 11. Naturally, the direct contribution to the amplitude is independent of Aq 
(up to the very slight dependence of the squark spectrum) , but since altering Aq can 
affect the chargino and neutralino spectra by altering the values of /i and B which 
give the correct electroweak breaking and also affect the running of the off-diagonal 
mass terms, it can have a significant impact on the indirect contributions. 

We now touch upon how sensitive these results are to the other parameters. Select- 
ing 114 < gives a lighter chargino, which in general increases the chargino dominated 
rate, but in practice we find that since to evade the experimental chargino bound we 
must then pick higher values for the soft masses mo and M1/2, this does not in fact 
lead to very different behaviour. 

Increasing tan P can both directly enhance the contributions from diagrams with 
a helicity flip on the chargino or neutralino line and also change the spectra in a 
complicated way directly and through /i^. In fact, the first of these effects seems 
most important, as seen in Figure 12, which is similar to Figure 9, but with tan/3 
set to 30. Here the chargino and neutralino contributions both increase by around 
a factor of three, as we would expect if the dominant contribution to these is via 
higgsino-neutralino mixed states. 

In conclusion of our discussion of the effects generated by a pair of A' couplings, we 
have found that the dominant contribution to the amplitude is that of the chargino 
followed closely by that of the neutralino, with the direct contribution from the squark 
mediated diagrams almost completely negligible. This clearly has important implica- 
tions in any attempt to bound R-parity violating couplings, since the only case where 
we can do a reasonably complete analytical analysis is irrelevant. In general, we can 
conclude that a choice of X[^jX2ij{MGUT) of order 10~^ is fairly safe for mo and M1/2 
of order lOOGeV, compared to a limit from the direct diagrams about two orders of 
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magnitude weaker. However, the various contributions are sufficently complicated 
that this can only really be regarded as an order of magnitude estimate, which may 
be evaded by appropriate choice of tan/3 and the various soft parameters. We em- 
phasise here that the amplitudes A are all directly proportional to the product of the 
two couplings, and as before it is straightforward to convert our figures showing the 
amplitude into bounds for a given point in parameter space. 

Having shown the general behaviour for the A' case, we now consider two alterna- 
tive scenarios. Firstly, we consider the case Xi3i{Mgut) = A23i(Mg(7t) = 0.01. Here 
we note that these couplings are sufficiently large to give a direct contribution which 
is close to the present experimental limit. 

We show the resulting contributions to the amplitude as a function of M1/2 in 
Figure 13, with tan/? = 10, tjiq = lOOGeV, Aq = 0, /i4 > 0, and Xi3i{Mgut) = 
X23i{Mgut) = 0.01. The apparently pathological behaviour of the direct contribution, 
which changes sign even though it does not depend on the stau masses and mixings, 
can be understood as follows. From equation ( ^.4|) , we see that the amplitude behaves 
as 

Mr ~ -iy ^ (5.21) 

ml 2m? ^ ^ 

For small M1/2, the slepton and sneutrino masses will be similar, and so the first term 

on the right hand side of this equation will be larger and the result positive; however 

as we increase M1/2, ttlI^ will grow larger than 2m|^ due to its greater electroweak 

radiative corrections in the RGBs, at which point the amplitude changes sign. We 

find this behaviour occurring when m,>^ and mg^ are around 230GeV and 165GeV 

respectively. Here, however, the chargino and neutralino contributions are still larger 

than that of the direct contribution, although much less so than in the case of the A' 

diagrams. 

In Figure 14 we show the same data as for Figure 13 except for tan (3 = 2, but with 
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Ml/2 held constant at lOOGeV, and varying ^40. This is an example of exact cancella- 
tion of the amplitude. Since we often find that the direct and indirect contributions 
to the amplitude are of similar magnitude and opposite sign, there are many cases 
where the cancellation can be very good, particularly if tan (5 is not too large (since 
if we increase tan j3 we increase the chargino and neutralino mediated contributions 
to the amplitude without greatly changing the slepton mediated contributions). 

Hence, our conclusion for the contribution to /i ^ 67 from setting A131A231 non- 
zero is ultimately that the new contributions from chargino and neutralino mediated 
diagrams dominate for large tan /3, but that there are some fairly subtle cancellations 
which become increasingly significant for smaller tan/3. Often the direct and indirect 
contributions have opposite sign, and hence it is not really possible to derive a reliable 
bound. 

The last scenario which we shall consider is that of Arl, which we generate by 
setting X\'2.i{Mqut) — y^i2i{MGUT) — 0.01. Again, we expect the contributions from 
the direct diagrams to be relatively close to the experimental limit. However, this 
case is rather different from the two previous ones because we no longer have a 
chargino contribution. Hence we would expect that the cancellations between direct 
and indirect effects might occur for larger tan/5. That this is the case is shown in 
Figure 15, in which we show the amplitude as a function of M1/2 with again tan {3 = 10, 
mo = lOOGeV, — {)^ jji,^> 0, and \i2i{Mgut) — ^i22{Mgut) — 0.01, showing that 
the cancellation ensures that the amplitude remains well below the experimental limit. 

We conclude this section by summarising our results. Firstly, we find rather 
different behaviour for the three different scenarios which we have considered. For the 
A' case, the chargino and neutralino mediated diagrams with fiavour violation through 
soft mass insertions dominate completely the direct contributions, giving very much 
tighter constraints, particularly for large tan /3. For the LR effects due to A couphngs 
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we find again that the chargino contribution dominates, but not overwhelmingly, and 
there can be large cancellations. For the RL case there are no chargino contributions, 
and the neutralino and direct effects are usually of comparable size and opposite 
sign. In all of these cases it is very hard to do much more than gain an idea of the 
likely size of the effects in question since the parameter space is complicated, but 
it is clear that the Ai3jA23i and especially X'njX^ij amplitudes will be much greater 
than those from the direct contributions alone, while the neutralino contributions 
to XijiXij2 are of the same order of magnitude as the direct effects. However, since 
there are so many possible cancellations between terms, it is essentially impossible 
to derive concrete bounds. The strongest reasonable statement is that, for the values 
we have considered for pairs of couplings at Mqut of AA ~ 10~^ and A' A' ~ 10~^ we 
expect contributions of order the experimental limit for a very light spectrum, with 
the branching ratio scahng as A^ or A'^. 

6 Conclusion 

In this paper we have given full RGEs for the MSSM with R-parity violation with 
the inclusion of all soft terms as well as all dimensionless couplings. We have given 
solutions to these equations for small R-parity violating couplings to relate values 
at the unification scale to those at low energy, and given the triviality bounds on 
couplings. 

The inclusion of R-parity violation in our superpotential through dimensionless 
terms allows the generation of lepton-Higgs mixing which leads to sneutrino VEVs 
and hence neutrino masses. We have presented a detailed discussion of sneutrino 
VEV generation, together with a calculation assuming universal soft terms respecting 
lepton number at the unification scale. This shows that the indirect generation of 
sneutrino VEVs through the running of the RGEs for the soft terms often leads to 
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larger effects than those derived directly from one loop diagrams, and hence that 
Ai33, A233, A^33 must all be small. Typically we find that values of Aj33 and A'j33 of 
order lO'^ and 10"^ at the GUT scale give masses to the corresponding neutrino of 
order hundreds to thousands of eV, although the exact value is quite dependent on 
the unification scale parameters. 

Similarly, R-parity violation at the unification scale can generate large lepton 
and quark flavour violation both directly through one loop diagrams and indirectly 
through the generation of off-diagonal soft masses and hence flavour violation in the 
scalar sector. We have studied the process fi 67, which we have shown to be very 
strongly affected by chargino and neutralino mediated diagrams. These typically 
dominate the direct contributions which had already been calculated, typically by 
several orders of magnitude for the case of the A' couplings, but there are often 
cancellations so that it is not possible to give precise bounds on the couplings from 
such processes. However, unless we invoke arbitrary cancellations, the typical size of 
such indirect effects on FCNC are likely to be the dominant constraint on the building 
of a model with non-zero R-parity violating couplings. 

Our main conclusion from both these calculations is that R-parity violation can 
generate large lepton flavour and lepton number violating effects through the running 
of the dimensionful RGBs, and that these effects are often much larger than those 
which are generated directly by the couplings themselves. However, it is virtually 
impossible to turn this statement into hard numerical bounds because of the number 
of free input parameters and the large cancellations between different contributions. 
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A Renormalisation Group Equations 

We now present our conventions and RG equations. The superpotential and soft 
potential are given by 



W = htQ3H2U3 + hbQsHids + KLsHies 
-'r-^KjkLiLjek + \'ijj.LiQjdk + 



Vsoit = ritQ3H2U3 + TjhQsHid^ + r/^LgifiCs + h.c. 

1 ^ 1 

~\"^CijkLiLjek + C^jj^LiQjdk + —C^jj^Uidjdk + h.c. 

+D4H1H2 + DlhH2 + h.c. (A.l) 

Here the notation is conventional, with the exception that we do not write the trilinear 
terms as Ax coupling since this is not convenient algebraically in the context of R- 
parity violation, and similarly we do not write the bilinear soft term as B^. It should 
be noted that we now label /i as so as to make clear the analogy between the R- 
parity violating bilinear terms Hi and n^. Indices indicate generation, and all gauge 
indices are suppressed. Raising and lowering generation indices has no particular 
significance. In our expressions we will use /i interchangably with /i4 and for 
r]^/h^, and we express the mass eigenstates as e*^"^-*, e'-^\ and similarly for quarks. 

The derivation of the RG equations is straightforward using standard techniques 



found for example in |2J, |3^. Here, we use the convention that all indices which 
do not appear on the left hand side of the equations are summed over for the three 
generations, and we define 
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where the sum is over all particles, and must include all appropriate degrees of free- 
dom factors from, for example, colour and generation. Weak hypercharge is here 



normalised so that Yq — 1/6. We use t — In//, where is the MS renormalisation 
scale. 



dt 



7 16 
+^rA„j33A^33 - hi,{-gl + ?,gl + —gl) 



dt 



dt 



,dX' 



+?>hkXmz^X!^^^-hr{?>gl + ?,gl) (A.3) 



IGtT _ XiyYiffi{XijyiYiXij}^ -\- XijYik^ljn ~l~ AjjYjA;,7^/j) 



' Imn 

+Aijfe/l^((^i3 + (^j3 + 25fc3 + 5jr^5kZ + (^iSf^fes) 
+ /lf,/lT-(3(5j3(5fc3A-33 — 35i35fc3A^-33) 

-Xijk{3gl + 35^2) 



167r^~^ — Aj^„(3A^^„Ajj-;^ + 2A^j-„A;^;^ + A-^;;.A;^- J 

''fXlmnXij^Ximn + '^^lmn^ijn\mk 

+Kjki^i3hl + ^j^hl + 5j3/i6 + 25fc3/l^ + 35j3(^fc3/i6) + SjsSkshbhrXiss 

-Kjki^gf + ^gl + Y^I) 

d\" 

1fi-7r2__!2^ — \" (\" \" _|_9\" \" _l_9\" \" ^ 
-■-O^ — '^lmn\\mn'^l3k '^\mk'^ljn "^^jn'^lmk) 



+A:;.,(25,3/i? + 25,3/^^ + 243/^6 



-A::,,(^^? + 8^3') (A.4) 
167r^^ = r7*(18/i? + /i^ + AU' + ALn') 



+ht{2hbritj + 2A„j3„C^3„ + 2X^^^C^^^j 



13 16 13 16 

+ 3^2 + yffi) + 2/it(-^i'Mi + 3^2m2 + —glM^) 
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r]b{l8hl + hl + + GA^gg^ + 2A^„3^ + X'^sJ + 2A'^„3^) 
+r]r{2hbhr + 2A^33A^33) + 2hbhtr]t + C'^ssi^^bX'^^^ + KX^nss) 

+4/lbA^„3C^„3 + 2/l;,A^3„C^3„ + 4/l;,A^„3C^„3 

7 16 7 16 

-Vbi^gl + ^92 + y^I) + '^hbi-glM, + ^glM^ + —glM^) 

rjr{l2hl + 3hl + 2A^33 + A^3„ + A^„3 + SAg^^^) 
+rib{6hT-hb + 6Ato33A^33) + Cmssi^hbX'j^^^ + hrXmsa) 

+2/l^Am3nCm3n + ^hrXmnsCmnS + ^h-rX'^^^C ^^^) 

-VriSgl + 3gl) + 2hr{3glM, + Sg^M^) (A.5) 

Cljk + Xij 
~l~2Ajj^C/^/j + 2XijkCirnn ~l~ 2Xijnf;Cijn) 

+Cijkhl{Si3 + Sj3 + 2Sk3 + 2Sj3Sk3 + 2Si3Sk3) 
+Xijkhr7]r{26i3 + 26j3 + 4^3 + ^js^feS + fe^fes) 
+hbrir(S5j35k3X[^S - S5i35k3X'j33) 
+2/i6/ir(3(^j3<^fc3C'-33 — ^Si35k3Cj^^) 

-Ci,k{3gl + 3gl) + 2X,jk{3glMi + 3glM2) 

A' fsA' r'' + \' c' + 2A' r'' 

lmn\ imn Ijk ' Ijn imk ' Imk ijn 

+^Kjn^lmk + '^Kmk^ljn + ^Kjk^imn) 

'^XiYnn{XifjinCijj^ -\- 2Xij^Cimn) 

'^Kmni'^Kmk^'ijn + ^Kjrfilmk) 

+CljkiSi3hl + 5j3h1 + 5j3hl + 26k3hl + Q5j3h^hl) 
+X'ij,,{26i3hrr]r + 26j3htr]t + 26j3hbr]b + 4:Sk3hbr)b + S^js^fesMft) 

+ ((^j343^rAi33?76 + 25j343^&^TCi33) 

+ + y^3') + 2A^,(^^?Mi + 3^,^M2 + ^^3'M3) 

'^ImnWmn^ljk ^"^Ijn^imk ^"^Imk^ijn 

'^^KjvPlmk + 4A^'^feC;'^„ + 2X!ljj^C'l^^ 

'^\mn{'^\mj^ink + '^\mk^ijn + ^Kjrfilmk + ^\nk^lmj) 

32 



+C'/j^i25i3h^,+25j3hl + 25k3hl) 

+\''jj,{A6i^htr]t + 4:6j3hbr]b + A6k?.hbr]b) 

-CU^gl + + KdlglM, + l&glM,) (A.6) 

'^2CmniCmnj + 25ij$, — SSijgfM^ 

'^'^Kmn^'kmn''^%Lk + '^^jmn'^kmrJ^\iLk 
+2Xii^Xjlnml^^^ + 2Ximl\jnlfn\^L^ 

+mlLj {Sishl + 5jshl) + 5i^5j^{2hlm]j^ + 2hlm%,^ + 2r)l) 

\i \i ™2 _|_ \i y 2 
^min'^mkn" ''QjQ^. '^mjn'^m.kn" ''QiQk 

'^'^^kim^kjrJ^Xndn + 2A^jj^.A^jj^.m^^j;,^ + 2C'^^JJ'^j^ 

+6is6jsi2h^tmjj^ + 2hlml^^^ + 2ril + 2/i2m2,^ + 2/i2m^3^3 + 2ril) 

'^imn'^kmn^UjUi^ ' '^jmn'^kmn'^UiUj, ' '^'^ikm'^jkn^dmdn ' imn^jmn 

-\5,,i - 25,,{^glMl + fglMi) 

'^^mni^mnk'^djdk '^^mnj'^mnk^didu 

'^'^^mni^mnk''^djdk + '^^Lnj^Lnk'^'didk 

''l'^\mi\nj''^Q,nQ„ + '^\nkiKikj''^'LrnL„ 
_L/1 \" \" 2 _L /I \" \" 

'^^'^kmi'^kni^dmdn + ^^rnki'^nkj''^' 



2 



' mni mnj ' mni mnj 
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2 4 Ifi 

- 25,,A-glMl + -^32M|) (A.7) 

In the above equations, we have not allowed for the inclusion of Yukawa terms 
other than those of the third generation, and we have also not presented dimensionful 
terms such as m\.H^ which can be generated by the mixing of Li and Hi. However, 
these may be trivially derived from the above by simply regarding Hi as L4, dropping 
the explicit hi, and h^. dependence of the equations, and setting 

-^343 = "-^433 = hr A433 = — /ife C343 = — C433 = 7/^ C433 = — ?7fe (A. 8) 

Adopting this notation, we find 

16vr2^ = ^i(3/i? - gl - 3^2) + 'iKki^^'jkiN + ^iki^kil^j 
16vr2^ = D,{?,hl-gl-?,gl) + ?,K^i\'^^iDj + Kui\,kiDi (A.9) 
+6htr]tfii + GC'i^iX'jkifJ'j + 2CikiXjkifJ'j + 2glMifii + 6glM2fJ,i 

We further note that although we have given the RGEs for the third generation 
Yukawas only, it is straightforward to use this notation to generate the full RGEs 
including general mass matrices and the full CKM dependence through 

AW = -hf, (A. 10) 

where hfj is the full Yukawa matrix for the down-type quarks, which can at low 
energies be found from the quark masses and the CKM matrix as defined, for example. 



in 1 22 1 . In general, the R-parity violating couplings will prevent all the mass matrices 



remaining diagonal even if they are initially so chosen, as discussed above. 

For completeness, we also give the RGE for mfj^jj^ , which is unchanged from that 
of the MSSM. 

-{2glMl + QglMl)+i (A.ll) 
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B Definition of Functions 



The expressions for the amphtudes in /i — > 67 presented earher employ a number of 



functions, most of which can be found in either ref. or in Appendix B of ref. ^8 
We repeat the definitions here for convenience. 

1 



F2{x) 

G{x) 
F{x) 
H{x) 
L{x) 



12(x - 1)4 

W-) 

X X 

1 

12(x - 1)4 
1 



[x — 6a; + 82; + 2 + Qx ln(a:;)] 



[2x^ + 3a;2 - 6x + 1 - Qx"^ ln(x)] 



\x 



2{x-lf 
1 

2(x-l)3 
Fi(x) +xF[{x 
1 



4a; + 3 + 2 ln(x)] 



[x^ - I- 2xln(x)] 

+ - 9x - 1 - 6x(x + 1) ln(x)] 



6(x- 1)5' 

F2(x) +a;F^(a;) 
1 

— rT[-17x^ + + 9x - 1 + 2^2/^ ^ j^^/^ni 

12(x — 1)^ 

F3(x)+xF^(x) 

^ -[x2 + 4x-5-2(2x+l)ln(x)] 



2(x - 1)4' 

^) 
1 



F4(x) + xF4(x) 



-Sx^ + 4x + 1 + 2x(x + 2) ln(x)] 



(B.l) 



2(x- 1)4' 

Here the prime simply indicates the first derivative. Certain asymptotic values of 
these functions are useful, as listed in Table 2 below. 



X 


Fi{x) 


F2{x) 


F,{x) 


Fi{x) 


G(x) 


F(x) 


H{x) 


L(x) 





1/6 


1/12 


— ln(x) 


1/2 


1/6 


1/12 


- ln(x) - 5/2 


1/2 


1 


1/24 


1/24 


1/3 


1/6 


1/60 


1/40 


1/12 


1/12 


00 


l/12x 


l/6x 


l/2x 


l/2x 


1/6x2 


ln(x)/2x2 


1/2x2 


ln(x)/x2 



Table 2. Values of the various functions in useful limits. 
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C Figure Captions 



Figure 1: Fcynman diagrams contributing to Lj — Hi mixing and iience to the 
generation of sneutrino VEVs. We do not show possible mass insertions on the hnes. 

Figure 2a: Absolute values of the sneutrino VEV (dashed lines) and corresponding 
neutrino mass (solid lines) as a function of m^., with parameters rrit — 175GeV, 
asiMz) = 0.12, tan/3 = 2 and 20, M1/2 = SOOGeV, Aq = 0, both signs of and 
Xi33{Mgut) = 0.01. 

Figure 2b: Absolute values of the sneutrino VEV (dashed lines) and corresponding 
neutrino mass (sohd hnes) as a function of m^,., with parameters rrit — 175GeV, 
a^^Mz) = 0.12, tan/3 =2 and 20, M1/2 = SOOGeV, Aq = 0, both signs of //, and 
A'i33(Mgc7t) = 0.001. 

Figure 3: Absolute values of the sneutrino VEV (dashed lines) and corresponding 
neutrino mass (solid lines) as a function of M1/2, with parameters rrit — ITSGeV, 
a^^Mz) = 0.12, tan/3 = 2, 20, mi, = SOOGeV, Aq = 0, both signs of and 
KssiMour) = 0.001. 

Figure 4: Absolute values of the sneutrino VEV (dashed lines) and corresponding 
neutrino mass (solid lines) as a function of Aq, with parameters rrit = ITSGeV, 
a3(Mz) = 0.12, tan/3 = 2 and 20, M1/2 = mi, = SOOGeV, both signs of /x, and 
A;33(Mgc7t) = 0.001. 

Figures 5: Diagrams contributing to A^^j^. The analagous contributions to have 
a left handed internal lepton line. All diagrams are presented without photon lines 
for clarity, and the crosses indicate helicity flips except where otherwise marked. 

Figure 6: Diagrams contributing to A^'^. There are no analagous diagrams for A^'^. 

Figure 7: Diagrams contributing to A'^^. The analagous contribution to A'^'^ is 
proportional to the electron mass in the chargino case and hence neglected, while for 

39 



the RL neutralino case right and left handed fields are interchanged. 

Figures 8: Examples of diagrams which contribute to A^j^ with the insertion of a 
sneutrino VEV and hence neutralino-neutrino or slepton-chargino mixing. 

Figure 9: Absolute values of amplitudes for — > 67 from direct R-parity viola- 
tion diagrams (dashed lines), neutralino mediated diagrams (dot-dashed lines), and 
chargino mediated diagrams (dotted lines) plotted against M1/2. We also show the 
total amplitude (solid line) and the experimental bound on the amplitude (horizontal 
solid line). Parameters are = 175GeV, as{Mz) — 0.12, tan/3 = 10, mo = lOOGeV, 
Ao = 0, /X4 > 0, and X[^^{Mgut) = A^n (A^gc/t) = 0.001. 

Figure 10: Exactly as for Figure 9, but varying mo and with M1/2 set to lOOGeV. 

Figure 11: Exactly as for Figure 9, but varying Aq with both mo and M1/2 set to 
lOOGeV. 

Figure 12: Exactly as for Figure 9, but with tan/3 = 30. 

Figure 13: Absolute values of amplitudes for /j, ^ from direct R-parity viola- 
tion diagrams (dashed lines), neutralino mediated diagrams (dot-dashed lines), and 
chargino mediated diagrams (dotted fines) plotted against M1/2. We also show the 
total amplitude (solid line) and the experimental bound on the amplitude (horizontal 
solid line). Parameters are rrit = 175GeV, a^^Mz) = 0.12, tan/3 = 10, mo = lOOGeV, 
Aq = 0, //4 > 0, and \i^i{Mgut) = A231 (Mgc/t) = 0.01. 

Figure 14: Similar to Figure 13, but varying Aq. Parameters are m^ = 175GeV, 
az{Mz) = 0.12, tan/3 = 2, M1/2 = mo = lOOGeV, Aq = 0, //4 > 0, and Xi^i{Mgut) = 

A23i(Mgc7t) = 0.01. 

Figure 15: Absolute values of amplitudes for — > 67 from direct R-parity violation 
diagrams (dashed lines), and neutralino mediated diagrams (dot-dashed lines) plot- 
ted against Mi/2- We also show the total amplitude (solid line) and the experimen- 
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tal bound on the amplitude (horizontal solid line). Parameters are rrit — 175GeV, 

asiMz) = 0.12, tan/3 = 10, mo = lOOGeV, Ao = 0, ^4 > 0, and Xi2i{Mgut) = 
Ai22(Mgc/t) = 0.01. 
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